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Abstract 


Investigation of annular obstructed azimuthally polarized Bessel-modulated Gaussian beam 
(QBG) in the focal region of high NA lens, based on vectorial diffraction theory. The numerical results 
show that the intensity distribution in focal region of the incident beam can be altered considerably by 
changing beam parameter (1) and introducing annular apodization (6). Beam parameter induces the focal 
splitting in transverse direction, while annular apodization leads to change in focal pattern along optical 
axis of the focusing system. More interesting, the focal splitting may be in continuous in certain case of 
incident beam propagating through aligned optical system which is suitable for application such as optical 
manipulation and optical trapping. 
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1. INTRODUCTION 


The laser beams with cylindrical symmetry in 
polarization have attracted very much for their 
interesting properties and applications called 
cylindrical vector beams. This kinds of beams can be 
generated by using active and passive methods (Zhan 
et al. 2009). Due to the polarization symmetry, the 
electric field at the focus has unique polarization 
properties. K.S. Youngworth and T.G. Brown 
calculated cylindrical-vector fields (Youngworth et al. 
2000), which shows that, in the particular case of a 
tightly focused radially polarized beam, in the focal 
region the polarization shows large inhomogeneities, 
while the azimuthally polarized beam has purely 
transverse field in focal region. Focus shaping 
technique is also reported by using generalized 
cylindrical vector beams (Gao et al. 2007), in which a 
generalized cylindrical vector beam can _ be 
decomposed into radially polarized and azimuthally 
polarized components. And a generalized cylindrical 
beam can be generated from a radially polarized or an 
azimuthally polarized light using a two-half-wave-plate 
polarization rotator. Recently, the tight focusing 
properties of a polarized beam and a circularly, radially, 
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azimuthally or linearly polarized vortex beam by a high 
numerical-aperture (NA) lens have been discussed. The 
study of coherent superposition of radially and 
azimuthally polarized components has also generated 
interest (Jia et al. 2005). The polarization propagation 
and light intensity distribution in focal region plays 
an important role in many optical systems. It was 
demonstrated both theoretically and experimentally. 
Optical intensity distributions in focal region have 
found applications in optical micro manipulation 
domain (Veerabagu et al. 2013). Pure phase-shifting 
apodizer effect in focal region is reported (Mei et al. 
2005). Nowadays, the Bessel-modulated Gaussian 
beams with quadratic radial dependence (QBG beam) 
have novel focusing properties in cylindrical coordinate 
system. Caron and Potvliege introduced QBG 
(Veerabagu et al. 2013) it has attracted much attention 
(Wang et al. 2008). The axisymmetric QBG beam are 
in zeroth-order, which is usually referred to as the 
axisymmetric QBG beam (Dorn et al. 2003), Young 
worth and Brown calculated cylindrical-vector fields, 
near the focal region of an aplanatic lens, and briefly 
discussed some applications (Youngworth et al. 2000). 
Zhan and Leger reported a focus shaping technique 
using generalized cylindrical vector beam. In our 
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knowledge, almost all QBG beams in previous papers 
are in scalar form without considering the polarization 
distribution of optical field. In this paper we 
numerically study the focusing properties of 
azimuthally polarized axisymmetric QBG beam with 
radial variance phase wavefront. Though the proposed 
method is simple and free from the complex 
interferometric methods used in the conventional 
techniques to shift the generated focal hole. The 
principle of the focusing radially polarized 
axisymmetric QBG beam with radial variance phase 
wavefront is given in Section 2. Section 3 shows the 
simulation results and discussions. The conclusions are 
summarized in Section 4. 


2. THEORY 


In the focusing system we investigated, 
focusing beam is azimuthally polarized axisymmetric 
QBG in radial [2, 9, 28, 32]. Therefore the cylindrical 
coordinate system (7,¢,0) the field distribution E(r,¢, 0) 
of the azimuthally polarized axisymmetric QBG beam 
at the plane z = 0 is written as, 


E,(r,9,0) = Ey (r,9,2 = 0)-[cos(¢(r))-n, | (1) 


Where ¢(r) is the polarization angle from azimuthal 
direction, and is the function of azimuthal coordinate 
for the azimuthally polarized axisymmetric QBG beam, 
were n; is azimuthal unit vector. 


2: 2 
E, (7,9,2=0)=J, [2 es9{ | (2) 


0 WD 


Where Jp denotes the Bessel function of order zero, wo 
is the waist width of the Gaussian beam, pt is a beam 
parameter which is complex valued in general. After 
some simple derivation [29], Eq. 2 can be rewritten as 


Bree-0)=1{ SP xf SA Noite) ©) 





Parameter W = @,/fo is called relative waist width, 
where r, is radius of incident optical aperture. NA is 
numerical aperture of the focusing system. @ is polar 
angle corresponding to radial coordinate. ¢(r) is the 
polarization characteristics of the focusing azimuthally 
polarized axisymmetric QBG beam and can be 
expressed as 


_tan (@)r (4) 


a ea tan (@ ) 


Where a = arcsin (NA/n) is the convergence angle of 
the focusing optical system. If the focusing optical 
system investigated is in air, i.e. refractive index n= 1. 
C is radial variance parameter indicating polarization 
radial degree the wave front phase distribution. Using 
the same analysis method as that in references [2, 30], 
the electric field in focal region of azimuthally 
polarized axisymmetric QBG beam is, 


E(r, 9.2) =e, +H,e,+£,é, (5) 
Elrsge)-24| [eal sin A G.I drsineptiteasidd (6) 


E(r,g2)=E(r,g2) 0 


Where A is a constant, J,(x) and J;(x) denote zero 
order and first order Bessel functions of the kind, k is 
the wave number k = 27/1, with i being the wavelength 
of illuminating beam. Were e;, e» and e, are unity 
vector. Substitute Eq. (3) and (4) into Eqs. (6), then 
substitute Eqs. (6) into Eq. (5), we can obtain the 
analytical total electric field in focal region. The 
optical intensity in focal region is proportional to the 
modulus square of Eq. (5), so the focusing properties of 
azimuthally polarized axisymmetric QBG beam can be 
investigated theoretically. 


3. RESULTS AND DISCUSSIONS 


Without loss of validity and generality, it is 
proposed that the numerical aperture of the focusing 
system NA = 0.95 and amplitude constant A = /. The 
focusing properties of azimuthally polarized 
axisymmetric QBG beam is investigated theoretically, 
for parameter C = 0 at different beam parameter (uw). 
Firstly, the intensity distributions in focal region of the 
azimuthally polarized QBG beam are calculated under 
condition of C = 0 and different u, and are illustrated 
in Fig. 1. In this paper, it should be noted that the unit 
of coordinates in all figures in this article is A~’. Focus 
refers to the maximum optical intensity peak, and focal 
shift denotes the movement of this focus. 


Fig. 1 shows the total intensity distribution of 
incident azimuthally polarized QBG in the focal region 
of high NA lens (NA = 0.95). It can be seen from Fig. 
l(a) in the focal region for w = 0, i.e. the system 
generates two focal spot with equal intensities in radial 
direction forms focal hole of focal depth (1.44) each. 
The focus (focal hole) in the focal region splits along 
optical axis on increasing parameter uw, namely focal 
split phenomenon occurs and is shown in Fig. 1(b) for 
ue = 4. On increasing uw = 6 continuously, the focal 
pattern changes considerably and illustrated in Fig. 
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1(c). There is two bright spot with two weaker spot on 
either side of bright spot in the focal region. We further 
increased uw = 9, in order to study the effect of beam 
parameter p and its evolution in the focal region is 
shown in fig. 1 (d). We observed that on increasing 


the generated focal hole splits along optical axis and 
forms four bright spot in the focal region, the bright 
spot becomes weaker and the weaker spot becomes 
stronger. 

















Fig. 1: Total intensity distribution of azimuthally polarized high NA lens (NA = 0.95) for (a) u=0, (b) u=4, 


(c) u=6 and (d) u=8 


Secondly, we analyzed numerically the effect 
of annular apodization for different parameter uw by 
keeping other parameter same as above. The effect of 
annular apodization on the incident radially polarized 
axisymmetric QBG beam is analyzed replacing E(@) by 
E(@)T(@) in Eq. (5). Where T(@) is given by 


0 for 0<@<0d-a 
(a=) F 
1 for 6:a<s@0<a 
(6) 
In order to study the effect of annular 
apodization we choose annular apodizer 6 = 0.85, by 


keeping other parameter same as above for different 
pupil beam parameter and it is shown in Fig 2. 





(a) n=0 (b) n=3 


play 
pia 





pia 
pid 





2(Aa) zy 





(c) u=4 (d) n=6 


pla 
pin 





() u=9 (h) 


pla) 
pia) 











Fig. 2: Total intensity distribution of azimuthally polarized high NA lens (NA = 0.95) with annular 








obstruction (6 = 0.85) (a) u=0, (b) u=3, (c) u=4, (d) u=6, (e) u=7, (f) u=8, (g) H=9 and (h) “=10 
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Fig. 2 (a) shows the total intensity distribution 
contour plot of incident annular obstructed azimuthally 
polarized axisymmetric QBG beam when i = 0, the 
focal region has only one focal hole (two focal spot in 
radial direction with FWHM (0.28A) and its focal depth 
of FWHM (4.06A) of each). Further increase in beam 
parameter p, results in focal splitting along axial 
direction at the focus and it is shown in Fig. 2 (b) for uw 
= 3, for u = 6, generates focal hole in the focal region, 
it is shown in Fig 2 (c). Fig 2 (d) for uw = 6, Fig. 2 (d) 
for uw = 7, Fig. 2 (d) for wu = 8, Fig. 2 (d) for uw = 9 and 
Fig. 2 (d) for uw = 10. It is observed from Fig. 2 on 
increasing , continuously, the focal pattern changes 
from one focal hole to two focal hole in the focal 
region. 


4. CONCLUSION 


Based on vector diffraction theory, tight 
focusing of a azimuthally polarized QBG beam by a 
high NA lens is studied theoretically. The total optical 
intensity distribution in the focal region is illustrated by 
numerical calculations without and with annular 
apodization for different beam parameter yw. The 
simulation result shows, on increasing beam parameter 
u the proposed system generates different focal pattern 
changes at the focus. We expect such a beam can be 
widely used in application such as optical trapping, data 
storage, biomedical imaging, laser drilling, and laser 
machining. 
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